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‘Advanced  concepts  and  various  approaches  on  the  transducer  compensation 
module  along  with  the  concepts  on  an  intelligent  flaw  detector  module 
are  presentedAln  this  report. 
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I.  SUMMARY 


Classification  Applications",  was  conducted  at  Drexel  University  during  the 
contract  year  1977.  This  part  of  the  research  program  serves  as  an  initial 
phase  in  the  development  of  a larger,  more  expanded  program  on  the  develop- 
ment of  detailed  specifications  for  an  intelligent  flaw  detection  system. 

■^An  intelligent  flaw  detection  system  is  a computer  controlled  ultrasonic 
inspection  system  capable  of: 

y indicating  the  correct  inspection  procedure  for  parts  being  inspected; 

} analyzing  the  transducer  and  system  operation  to  insure  that  meaningful 
data  is  acquired;  ar~A 

} interpreting  the  data  and  compensating  for  any  system  variations  that 
may  occur  due  to  changes  or  aging  of  the  transducer  and/or  couplant 
variations  or  any  other  inconsistencies  in  the  system.  ^ 

The  problem  of  transducer  compensation  was  the  principal  subject  of  study 
during  this  initial  phase.  Phase  I.  The  results  obtained  using  a deconvolution 
algorithm  were  not  as  accurate  as  desired  due  to  both  "noise"  and  nonlinearities 
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II.  PREFACE 

A.  A feasibility  and  design  study  on  "An  Adaptive  Automatic  Accept- 
Reject  Module  for  Transducer  Evaluation  and  Potential  Flaw  Classification 
Applications",  was  conducted  at  Drexel  University  under  U.S.  Navy  Contract  No. 
N68335-77-C-0556.  It  was  found  that  linear  mathematical  modeling  used  for 
describing  and  evaluating  various  transducer  transfer  functions  was  useful 
for  solving  the  very  complex  problem  of  transducer  compensation,  but  in  a very 
limited  way.  A "GO/ NO-GO"  acceptance  criterion  for  the  transducers  had  to  be 
specified  before  employing  compensation  mathematical  modeling  analysis. 

B.  This  concept  made  use  of  a list  of  specific  transducer  performance 
features  as  measured  from  the  ultrasonic  signal  features.  The  feature  list 
consisted  of  such  items  as  transducer  center  frequency,  6 dB  down  frequency 
bandwidth,  near  field  value,  dB  variation  within  near  field,  a measure  of  beam 
symmetry,  etc.  The  "GO/NO-GO"  concepts  compared  features  with  a list  found 
suitable  from  advanced  flaw  or  material  analysis.  Percentage  deviations  of 
Individual  items  on  the  list  are  determined  from  feedback  mechanisms  associated 
with  the  classification  algorithm  development. 

C.  This  part  of  the  program  serves  as  an  initial  phase  in  the  development 
of  a larger,  more  expanded  program  on  the  development  of  an  intelligent  flaw 
detection  system,  the  transducer  playing  a tremendously  critical  role  in  the 
development  of  the  complete  system.  Basically,  an  intelligent  flaw  detection 
system  is  a computer  controlled  ultrasonic  Inspection  system  capable  of  in- 
dicating correct  inspection  system  procedures,  specifying  and  checking  trans- 
ducer and  couplant  performance,  obtaining  useful  ultrasonic  data,  providing 
compensation  for  any  system  and  transducer  variation  and  finally,  employing 

a classification  algorithm  for  determining  flaw  size  and/or  flaw  character- 
istics and/or  material  variations.  The  system  will  operate  on  an  interactive 
basis,  allowing  the  operator  to  communicate  with  the  computer  in  conducting 
the  ultrasonic  examination. 
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VI.  INTELLIGENT  FLAW  DETECTION  SYSTEM 


A.  INTRODUCTION  AND  MOTIVATION. 

1.  The  function  of  an  Intelligent  flaw  detector  is  to  combine  the  speed 
and  accuracy  of  a minicomputer  with  the  data  acquisition  capability  of  a so- 
phisticated ultrasonic  inspection  system  and  various  advanced  signal  processing 
techniques  into  a single  system  capable  of  conducting  an  inspection  sequence 
including  data  evaluation  with  a minimum  amount  of  human  intervention  and 
system  variation.  Total  intelligent  flaw  detection  systems  are  necessary  and 
needed  because  of  inherent  limitations  of  conventional  ultrasonic  systems.  For 
example,  the  potential  success  of  conventional  ultrasonic  inspection  techniques, 
to  a large  extent,  is  operator  dependent.  An  ultrasonic  inspector  while  per- 
forming an  inspection  is  usually  required  to 

a.  Select  the  proper  frequency,  type  and  size  of  transducer  from  a 
wide  variety  of  transducers, 

b.  Select  the  most  appropriate  inspection  technique, 

c.  Carry  out  the  inspection, 

d.  Evaluate  the  test  results. 

2.  Evidently,  the  success  of  such  an  ultrasonic  inspection  critically 
depends  on  the  experience,  skills,  working  conditions  and  psychological  factors 
of  the  operator.  Many  times  different  operators  select  different  inspection 
techniques  for  inspecting  the  same  part  and  come  up  with  contradictory  results. 
It  ia  a rather  well  established  fact  that,  compared  with  a computer,  humans  are 
quite  poor  at  handling  and  evaluating  the  ultrasonic  test  data.  In  other  words, 
even  a skilled  operator  cannot  positively  assess  flaw  criticality  and  its  quan- 
titative relationship  with  fracture  mechanics  theories,  etc. 

3.  Another  problem  with  ultrasonic  inspection  today,  is  associated  with, 

transducer  performance  evaluation  and  assuring  the  repeatability  of  the  ultra- 
sonic test.  It  would  be  desirable  to  compare  output  on  a given  day  with  a 
particular  instrument  and  operator  with  results  obtained  from  some  earlier  data 
acquisition  period,  A repeatable  ultrasonic  teat  result  is  definitely  important 
and  needed  in  order  to  solve  many  complex  problems  of  flaw  character! ration  and 
classification  using  pattern  recognition  algorithms.  Preliminary  research 
indicated  that  an  Invariable  ultrasonic  input  function  from  a broad  band  trans- 
ducer is  highly  desirable  to  obtain  the  highest,  reproducible  and  meaningful 
index  of  performance  for  flaw  classification  purposes  employing  sophisticated 
pattern  recognition  algorithm.  This  generally  is  not  possible  due  to  various- 
reasons:  for  example,  changing  the  controls  on  pulser-receiver  units,  trans- 

ducer aging,  transducer  wear  and  tear,  couplant  thickness,  pressure  applied 

to  hold  transducer  during  contact  testing,  etc. 

. The  problems  outlined  above  are  particularly  addressed  in  our  proposal 
on  "An  Intelligent  Flaw  Detector  Specification  and  Feasibility  Study  for  the 
I'  . Navy",  with  an  accent  on  automat  ion  utilizing  microprocessor  technology. 

»t  ures  of  an  intelligent  flaw  detector  include: 

. Increase  of  testing  speed  (reduced  downtime) 

. Klimination  of  human  failures 
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Automatic  transducer  evaluation 
Improvement  of  test  reproducibility 

Establishment  of  optimum  test  conditions  with  respect  to  flaw  detection 
end  eva lust ion 

Comprehensive  recording  of  results 
Automatic  evaluation  of  test  results 

Merger  of  ultrasonic  technology  with  microprocessor  technology 

Capability  of  carrying  out  advanced  analysis  in  pattern  recognition 
and  computer  learning  analysis 

5.  The  Intelligent  flaw  detection  system  will,  therefore,  be  a computer 
controlled  ultrasonic  Inspection  system  capable  of  Indicating  the  correct 
inspection  procedure  for  the  part  being  inspected,  analyzing  the  transducer 
and  operation  of  the  system  to  insure  that  meaningful  data  Is  being  acquired, 
interpreting  the  data  and  compensating  for  any  system  variations  that  may 
occur  due  to  change  or  aging  of  transducer,  change  of  couplant  or  any  other  change 
in  the  system.  The  problem  of  transducer  compensation  was  the  subject  of 
study  during  Phase  I (contract  year  1976-77)  and  presented  next. 

B.  TRANSDUCER  COMPENSATION  MODULE 

1.  An  investigation  on  the  subject  of  transducer  performance  evaluation 
entitled,  "An  Adaptive  Accept-Reject  Module  for  Transducer  Evaluation  and 
Potential  Flaw  Classification  Applications",  was  carried  out  at  Drexel 
University  under  Contract  Number  N68335-77-C-0556  during  the  contract 

year  1976-77.  The  motivation  for  carrying  out  this  work  task  was  guided  by 
one  of  the  final  goals  in  the  Intelligent  Flaw  Detection  Scheme,  namely,  that 
of  flaw  classification  with  minimum  human  intervention.  Research,  to-date 
in  the  area  of  flaw  classification  using  pattern  recognition  algorithms 
indicates  that  it  is  very  essential  to  use  similar  ultrasonic  input  functions 
to  obtain  meaningful  results.  It  may  be  pointed  out  that  a number  of  variables 
in  an  ultrasonic  system  such  as  changing  the  controls  on  pulser-receiver  units, 
transducer  aging,  transducer  wear  and  tear,  couplant  thickness,  etc.,  all 
affect  the  RF  waveform.  A transducer  compensation  module  was,  therefore, 
envisioned  which  could  eliminate  such  critical  parameters  as  type  of  trans- 
ducer, couplant  thickness,  aging  of  transducer,  day-to-day  variations  in 
pulser  etc. 

2.  A portion  of  the  research  carried  out  during  the  initial  phase  was 
presented  in  the  form  of  a research  paper,  entitled,  "On  Selected  Aspects  of 
Phase  Analysis  in  Ultrasonic  Testing",  which  was  presented  at  the  11th  Symposium 
on  Nondestructive  Evaluation,  April  20-22,  1977  at  San  Antonio,  Texas  and  is 
included  in  Appendix  A for  reference  purposes.  Mathematical  details  on  the 
subject  of  transducer  compensation  using  a deconvolution  algorithm  and  the 
concept  of  "OO/NO-OO"  acceptance  criterion  are  explained  in  the  above  paper. 

M.'  t hemat  leal  details  were  implemented  on  a PDP-11/05  minicomputer  using  a 
111  emotion  fllOO  analog  to  digital  converter,  Tektronix  4015  graphics  display 
package,  Aerotech  UTA-2  ixtlser  and  gate  and  Tektronix  7704  oscilloscope.  The 
mathematical  details  of  the  deconvolution  algorithm,  although  fairly  simple 
in  concept,  does  cause  various  implementation  problems  on  a computer.  One 
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has  to  surmount  such  problems  as  division  by  zero,  taking  care  of  phase  angles 
and  complex  division  etc.  The  computer  program  vas  checked  against  theoretical 
transfer  functions.*  Excellent  results  vere  obtained,  strongly  supporting  the 
possibility  of  elimination  of  ultrasonic  Bystem  variables  as  mentioned  previously. 
The  usefulness  of  deconvolution  algorithm  using  real  flaws,  however,  w»b  very 
much  limited.  The  problem  turned  out  to  be  much  more  complicated  using  con- 
ventional ultrasonic  equipment  and  procedures  due  to  the  following  reasons. 

a.  Quite  often,  the  pulser-receiver  characteristics  are  dependent  upon 
each  other.  The  RF  waveform  changes  when  attenuation  settings  on  the  pulser  are 
changed.  It  was  found  that  the  pulser  needed  electronic  modifications  in  order 
to  separate  its  characteristics  from  the  receiver.  Technical  discussions  in 
that  regard  are  under  way  with  several  instrument  manufacturers.  Separation  is 
crucial  in  establishing  reliable  compensation  procedures  and  meaningful  transfer 
functions  of  various  components  in  an  ultrasonic  test  system. 

b.  The  deconvolution  algorithm  has  the  underlying  assumption  that  the 
ultrasonic  system,  in  this  case,  is  a linear  system.  Various  components  of  the 
ultrasonic  system  are  shown  in  Fig.  3 of  Appendix  A.  The  discrepancy  in  the 
results  obtained,  therefore,  can  be  attributed  to  the  nonlinearities  present  in 
the  system.  It  was  then  decided  to  establish  a range,  if  any,  in  which  an 
ultrasonic  system  acts  as  a linear  system  and  develop  a "GO/NO-GO"  acceptance 
criterion  on  certain  signal  features,  as  an  example,  center  frequency,  6dB 
down  bandwidth,  pulse  duration,  number  of  cycles,  etc.  along  with  the  mathe- 
matical deconvolution  algorithm  mentioned  previously.  It  was  found,  however, 
that  implementation  of  such  a program  was  very  tedious  and  could  not  be  carried 
out  due  to  insufficient  feedback  from  the  algorithm  and  the  large  number  of 
variables  involved.  Real  test  flaws  are  required  to  make  the  "GO/NO-GO" 
acceptance  criterion  a reality. 

3.  Although  somewhat  discouraged  by  this  approach,  it  was  decided  to  re- 
consider the  overall  goals  of  the  intelligent  flaw  detector  program.  The  trans- 
ducer compensation  concept  could  be  further  utilized  once  the  overall  specifics- 
tions  for  an  intelligent  flaw  detector  were  determined.  Gome  of  the  important 
design  considerations  for  such  a system  are  reviewed  below. 

1* . First,  a general  purpose  program  which  would  guide  an  operator  through 
selecting  a proper  probe,  performing  the  inspection  on  any  part  and  evaluating 
the  test  result  was  evaluated  for  its  potential  application  in  our  project.  A 
careful  engineering  analysis  revealed  that,  such  a computer  program  would  be 
rather  tedious  and  unsuitable  for  solving  simple  inspection  problems,  besides 
the  implementation  into  a microprocessor  as  a field  instrument  may  become  very 
difficult,  if  not  impossible.  Keeping  in  mind  the  final  goal  of  a microprocessor 
inspection  module  for  every  different  problem,  a more  basic  approach  was  adopted. 
Flowchart  concepts  in  Figure  1*  present  an  initial  design  for  an  intelligent 
flaw  detection  system  with  emphasis  on  the  first  step  of  transducer  selection. 

The  cumulative  ultrasonic  inspection  problem  can  be  divided  into  several  cate- 
gories, e.g.,  thickness  measurement,  flaw  detection,  flaw  classification  in 


’•'For  ions  of  this  work  will  be  presented  at  the  Spring  meeting  of  ASNT  in 
' w Orleans  in  n paper  by  J.  Rose  and  M.  Avioli  on  "Transducer  Compensation  of 
Val.i’  in  Flaw  Classifications"  and  is  attached  for  reference  purposes  in  Appen- 
dix R. 

* No  t e that  the  module  concept  flowcharts  can  be  assembled  into  one  large 
flowchart  by  matching  numbers  contained  on  the  individual  pages. 
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metals,  bond  Inspection,  composite  material  inspection,  etc.  Each  of  the 
above  mentioned  categories  are  treated  as  different  problems  as  shown  in  the 
charts  and  may  result  in  a different  microprocessor  module.  These  flowcharts 
are  based  on  general  ultrasonic  concepts  and  can  be  modified  for  each  part  to 
be  inspected.  A flaw  detection  microprocessor  module,  for  example,  may  be 
designed  for  use  on  several  different  parts.  As  mentioned  previously,  the 
flowcharts  presented  here  are  preliminary  and  rather  general  in  approach  and 
are  meant  only  to  explain  the  concept  of  an  intelligent  flaw  detector  in  prin- 
ciple. Three  problems,  namely  that  of  thickness  measurement,  flaw  detection 
and  flaw  classification  in  metals  are  presented  in  order  to  explain  the  prin- 
ciple. Similar  diagrams  could  be  made  for  solving  various  other  important  NDE 
problems.  Flowcharts  are  presented  for  thickness  measurement,  flaw  detection, 
and  flaw  classification  problems.  The  intelligent  flaw  detection  system  will 
operate  on  sn  interactive  basis  allowing  the  operator  to  communicate  with  the 
computer  in  conducting  the  ultrasonic  examination.  The  operator  in  this 
inspection  procedure  has  to  determine  the  inspection  mode,  determine  the 
material  and  thickness  class,  etc.,  and  physically  carry  out  the  inspection. 
Notice  that  the  algorithm  will  aid  the  operator  in  appropriate  transducer 
selection  and  performance  evaluation.  In  the  case  of  the  flaw  detection  and 
thickness  measurement  problems,  transducer  compensation  is  neither  critical 
nor  required  hence  simplifying  the  examination  procedure.  Moreover,  the  reject 
transducers  from  flaw  classification  could  probably  be  used  for  thickness 
measurement  producing  possibly  a considerable  economic  savings.  Hence,  the 
overall  problem  of  transducer  compensation  and  acceptance  could  be  considerably 
simplified.  This  will  not  only  result  in  the  simpler  algorithm  and  hence  more 
economical  microprocessor  modules,  but  also  considerable  savings  in  inspection 
time.  It  may  also  be  pointed  out  that  automatic  couplant  selection,  couplant 
check  and  data  evaluation  will  be  possible  using  an  intelligent  flaw  detection 
algorithm.  It  is  worthwhile  mentioning  that  many  of  the  boxes  in  the  accom- 
panying diagrams  denoting  'D',  will  be  determined  in  Phase  II  and  Phase  III 
of  the  research  program. 

5.  In  the  flowchart  describing  the  flaw  classification  program,  the 
transducer  compensation  problem  is  the  most  important.  The  following  is  pro- 
posed and  being  evaluated  at  the  present  time  in  order  to  solve  the  very  com- 
plex problem  of  transducer  compensation.  Out  of  a feature  vector  containing 
several  Important  features  describing  the  transducer,  only  relevant  ones  will 
he  evaluated  pertaining  to  the  problem  along  with  the  deconvolution  algorithm. 

o.  Table  I illustrates  a feature  vector  that  may  be  used  for  transducer 
characterisation  purposes  when  solving  a more  general  kind  of  inspection  prob- 
!‘m.  For  solving  a specific  problem,  however,  only  a few  relevant  features 
may  be  evaluated.  For  example,  for  the  thickness  measurement  problem  the 
Tv-atures  may  be  pulse  duration  and  center  frequency.  The  specifications  on 
these  features,  for  this  particular  problem,  may  be  within,  say,  10-15  per- 
cent. For  flaw  detection  problems  such  features  as  pulse  duration,  center 
frequency,  Fresnel  zone  and  beam  angle  of  divergence,  could  be  evaluated,  the 
specs  again  are  not  very  critical.  For  a problem  like  23  flaw  classification* 
study , however,  the  feature  vector  may  contain  6db  down  frequency  profile 
bandwidth,  number  of  depressions  in  the  U-6  MHz  region,  kurtosis,  number  of 


•dose,  J.L.,  "A  23  Flaw  Sorting  Study  in  Ultrasonics  and  Pattern  Recog- 
nition," Materials  Evaluation,  July  1977. 


A 


i—il  1 1*0 


NAEC-GSED-121 


i: 


it 

I 


l 

1. 


f 


| 


I f 


Page  10 

peaks , etc.,  besides  the  features  nentioned  above.  The  specs  now  become  very 
critical  for  this  problem  and  may  be  within  5 percent  if  only  a rigid  "GO/NO-GO" 
transducer  acceptance  criterion  were  to  be  used.  The  spec  limits  on  these 
features  may  be  related  by  using  a feedback  mechanism  from  the  deconvolution 
algorithm,  in  order  to  establish  a linear  range  for  the  algorithm  to  work. 
Further  work  in  this  area  will  be  carried  out  during  PhaBe  II  of  the  Intelligent 
Flaw  Detection  Module  Program. 
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Table  I - A Typical  Feature  Vector  for  Transducer  Condensation  Module 


Flaw 

Detection 


Thickness 

measure- 

ment 


Flaw 

Sorting 

Problem 


Pulse  Duration 
Center  Frequency 
Fresnel  Zone  Length 
Beam  Angle  of  Divergence 
6 db  Down  Center  Frequency 
Number  of  Peaks  in  the  Waveform 
No.  of  Spectral  Depressions 
6 db  Down  Frequency  Profile  Bandwidth 
Amplitude  Time  of  Pulse 
Rise  Time  of  Punse 


Bond 

Strength 

Problem 


Fall  Time  of  Pulse 
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VII.  CONCLUSIONS 


In  conclusion , the  transducer  compensation  problem  has  been  found  to  be 
much  more  difficult  than  originally  anticipated  if  a generalized  solution 
technique  were  desired.  The  problem  will  be  approached  from  the  viewpoint 
of  a total  intelligent  flaw  detector  program  with  compensation  treated  in  the 
future  in  a specific  manner  depending  on  the  problem  being  studied.  This 
report  presents  a first  step  towards  such  an  intelligent  machine.  Flowcharts 
are  necessary . These  charts  could  be  used  as  a teaching  aid  for  training 
inspectors,  for  developing  laboratory  facilities,  for  establishing  a specific 
solution  for  many  Navy  problems.  Automatic  recording  and  data  interpretation 
are  feasible  with  such  a system.  Flaw  classification  or  further  analysis  of 
any  kind  that  may  be  deemed  necessary  can  be  carried  out  at  a later  date. 
Transducer  characterization  is  very  problem  oriented,  resulting  in  consider- 
able savings  in  time  and  money  if  handled  carefully.  Couplant  checks  are  made 
automatically  from  time  to  time.  Algorithms  will  be  designed  for  minimum 
human  intervention,  proper  transducer  and  inspection  technique  selection, 
data  acquisition,  transducer  compensation  and  data  evaluation.  Although 
these  algorithms  will  be  developed  on  PDP-11  minicomputer  and  Biomat ion  8100 
analog-to-digital  converter,  an  initial  study  during  Phase  III  of  the  program 
undertaken  for  implementation  into  microprocessor  technology. 


ACKNOWLEDGMENTS  - Thanks  are  given  to  Krautkramer  Branson , Inc . , for  their 
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TRIGGER  OFF  THE  FIRS' 
ECHO  TO  AVOID  WATER 
RATH  AMD  CALIBRATION 
PROBLEMS 


pluc  i:j  thickness 

RANGE  AND  MATERIAL 
FDR  TRANSDUCER  SELECTION 


SCAN  THE  COMPONENT  AS 
REQUIRED  BY  CODE  - D 


RECORD  ALL  THE  DATA,  POSITION 
AMPLITUDE,  LOCATION,  ETC. 
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Article  1.  "ON  SELECTED  ASPECTS  OF  PHASE  ANALYSIS  IN  ULTRASONIC  TESTING" 
Joseph  L.  Rose  and  Gurvlnder  Pal  Singh,  presented  at  the  11th 
SjraposluB  on  Nondestructive  Evaluation,  April  20-22,  1977  at 
San  Antonio,  Texas. 
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ON  SELECTED  ASPECTS  OP  PUASE  ANALYSIS  IN  ULTRASONIC  TESTING 


Jaeeph  L.  Rot*  and  Curvlnder  Pal  Stack 
Draxal  Unlvaralty 
Philadelphia,  Pennsylvania 


taladtl  aapacta  of  Fourier  tfinifoi*  phaa*  analysis  associated  with 
ultrasonic  Inspection  are  reviews!  In  thta  paper.  Phase  signature  tech- 
» alC«tes  that  nay  be  useful  la  flaw  characterisation  and  classification  work 

cowered.  Hatheaatical  details  of  the  deconvolution  processes  are 
•*••*•*•*  wl,h  eaphasts  on  transducer  compensation  software  and  extensions 
(o  psttern  recognition  and  learnlni  machine  analysis.  In  particular,  four 
topics  are  discussed:  ahlftlni  characteristics,  deconvolution,  phase 
signature  techniques,  and  value  of  phase  in  learning  machine  analysis.  In 
•"■■•tyi  fundamental  concepts  in  phase  analysis  are  studied  that  will  en- 
ehle  researchers  to  advance  the  state-of-the-art  In  ultrasonic  inspection, 
•MN**  analysts,  and  lostrunantatloa  design  and  selection. 


P 


l 
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A.  INTRODUCTION 

Until  recently,  ultrasonics  has  been  used  as  an 
effsctlve  nondestructive  flaw  detection  tool. 

It  could  be  used  quite  readily  to  locate  flaws 
la  heeogeeeoua  Isotropic  ascertain.  Rut  with 
the  advent  of  composite  notarial*,  bonded 
materials,  the  tremendous  Increase  In  the 
conplenlty  of  structural  designs,  and  the 
conservation  of  naterlal  and  fuel,  mere  flaw 
detection  Is  no  longer  sufficient.  In  fact, 
it  Is  becoming  Increasingly  difficult  to  de- 
fine "flaw".  This  Is  also  contributed  to  the 
electronics  and  Instrumentation  technological 
advances.  The  question  now  becomes  "la  this 
flaw  detrirental  to  the  effectiveness  of  the 
structure?"  Before  this  question  can  be 
•nswered,  however,  it  la  necessary  to  deter- 
alne  as  many  characteristics  of  the  flaw  as 
possible  for  failure  analysis. 

"Ultrasonic  Spectroscopy",  a nare  given  by 
®*,l(<'e  [II,  has  been  used  to  obtain  more 
Information  about  flows.  Additional  details  on 
spectral  analysis  can  be  found  in  (f,J,4|.  It 
hee  also  been  used  to  study  transducer  wave— 
fntms,  measure  thickness  of  thin  materials, 
measure  bond  strength,  inspect  co-i-oslte 
ascertain  and  to  study  attenuation  and  micro- 
structure  in  a material.  Work  to  date  has 
Tailed  quite  heavily  on  amplitude  versus 
frequency  analysis,  with  little  attention 
pis  ed  i'll  the  values  of  phase  analysis.  The 
put  pose  of  this  paper  Is  to  outline  selected 


aspects  of  phase  analysis  that  might  be  useful 
la  advancing  the  state-of-the-art  In  ultnjo:ic 
tasting.  Tour  selected  topics  In  phase  analysis 
are,  therefore.  Included  below. 

Before  solving  a particular  problem.  It  is 
strtiury  to  understand  the  phase  profile 
characteristics. 

f(t)  * tine  domain  representation  of 
an  ultrasonic  signal 

Kw)  frequency  of  an  ultrasonic  signal 

Kw)  • |r<w)|  m R(u)  y jl(«) 

♦(w)  • tan*1  I(w)  ? 

I R(u>  J 

#(w)  •*  Phase  Angle 

B.  DISCUSSION 

1.  Time  Shifting  Characteristics 

The  subject  of  the  time  shifting  characteristic 
associated  with  the  pulse  locution  within  an 
electronic  gate  and  Its  affect  on  phase  ancle 
la  reviewed  in  this  section.  Although  details 
of  the  procedure  are  straight  forward  from  a 
mathematical  point  of  view,  they  ate  Included 
here  for  reference  purposes  to  remind  that  t‘e 
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characteristic  la  laportaac  and  to  provide  a 
guideline  tor  coping  with  this  shift. 

Lot  f(t)  represent  the  ultrasonic  pulse  shape  . 

la  the  tlae  doaaln  and  F(w)  Its  Fourier  trans- 
ton.  Let  f(t)  he  shifted  by  a tlae  constant ( so 
that  f(t-C)  * l(t).  The  function  g(t)  trans- 
forms to  C(w)  • F(u)  s”'3“Q  in  accordance  with 
the  Fourier  transfon  property.  Again, 

6(a)  - F(u)'1>K 

- F(a)[coo(wc)-3  Sla(at)] 

|C(«)  | - 0 F («)/*{ Cos1  (wg) )+| F(a) l 1 ( Sin1  UO  ^ 

- | F(w)  | [Coa2(aO>Sln1(uC)]‘‘ 

- | P(w) | (1) 

which  deaoostracss  that  the  pulse  shift  does 
not  alter  the  aagnltuda  of  the  Fourier  transfora. 
Tlae  shifting,  however,  results  In  a change  In 
the  phase  angle  as  shown  below 

F(a)  - | F(a) | ,!♦<“> 
whereas  tlae  shifted  pulse  reeults  la 

C(u)  - F(u)  e"J,K 


•r 

6(a)  - | F(w)  | 

- ) F(a) ) «) 


6(a)  • |c(w)  | a^“*  (3) 

Coablnlng  (2)  and  (3),  we  get 
1 6(a) | eJt(->-|F(-)|  <) 

froa  equation  (1) 

|c(a) | - |F(w)| 

♦ (a)  “4(h)  - w(l2k* 

v*>*rt  k *•  *°  Integer  and  Introduced  due  to 
the  periodicity  of  the  function. 

Hence  the  new  phase  angle  due  to  pulse  shift 
la  given  by  the  old  phase  angle  minus  mQ  if 
the  pulse  1,  moved  to  right  (away  from  origin). 
Similarly,  If  the  pulse  is  soved  to  the  left- 
hand  side  inside  the  gate  (cowards  the  origin), 
sew  phase  angle  la  given  by 


d(a)  • ♦(„)  a i 2kr 

In  other  words, 

♦new  " ♦old  * 1 2k* 

i app  u.  rlelely 

' * - 4 , . v k-> 

new  old 

>iKn  depends  upon  the  new  location  of 
:hr  puls*  as  compared  cr  the  'Id  one. 


A coaputer  ptograa  using  the  fast  Fourier  trans- 
fan  was  used  to  analyse  the  gated  pulse  in  the 
tlae  doaaln  to  obtain  the  Fourier  gpeccrua  and 
phase  angle  va.  frequency  data.  Variations  of 
phase  angles  due  to  roveaent  of  the  pulse  with- 
in the  gate  were  studied.  It  was  observed  that 
two  conditions  Bust  be  satisfied  in  order  to 
attain  the  above  results.  They  are,  namely, 

(1)  band  width  » 2 x center  frequency, 

aad 

(11)  that  the  shifted  (new)  pulse  lies  on 
the  sane  side  of  &■*— * as  the  old  one.  Satis- 
faction of  the  above  two  conditions  arises, 
naturally,  if  one  coneldcrs  the  discrete  Fourier 
transform  (5).  Saaple  results  are  illustrated 
la  Fig.  1. 


It  night  be  pointed  out  that  phase  profiles 
that  are  nonlinear  could  come  about  as  a 
result  of  examining  the  results  froa  multiple 
flaw  situations.  Additional  work  is  current- 
ly being  undertaken  on  this  project  to  com- 
bine both  theoretical  aod  experimental  aspects. 

2.  Deconvolut'"  n 

The  subject  of  transducer  compensation  is 
treaandously  critical  in  many  problems  of  ultra- 
sonic Inspection.  It  would  be  desirable  to 
coapare  output  on  a given  day  with  a particular 
instrument  and  operator  with  results  obtained 
froa  soae  earlier  data  acquisition  period.  A 
basic  deconvolution  model  la  outlined  in  thla 
papar  that  contains  some  hope  for  solving  this 
difficult  transducer  compensation  problem.  It 
Wight  be  pointed  out  that  results  to  date  do  not 
provide  us  with  a complete  solution  to  the  pro- 
blem because  of  the  nonlinearities  Inherent  in 
an  ultrasonic  system  as  well  as  the  available 
ultrasonic  input  pulses  that  must  be  used.  The 
future,  however,  shows  promise  for  a transducer 
compensation  algorithm  and  deconvolution  pro- 
vided this  mathematical  process  is  coupled  with 
an  acceptance  criteria  of  an  ultrasonic  trans- 
ducer. A "go"  "no  go"  acceptance  criteria  on 
certain  signal  features,  as  an  example,  center 
frequency,  fca3  down  band  width,  pulse  duration, 
and  so  on,  must  be  developed. 

A simple  deconvolution  model  is  presented  in 
Fig.  2.  For  a linear  system 


f(t)  * g(t)  - h(t) 
where  * indicates  convoluted  with 
or 


h(t)  - j f(t)  g(t-C)  d( 
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TIm  evaluation  of  this  Integral  caa  ba  carried 
Mt  In  a computationally  efficient  fashion  in 
frequency  doaala  atnca  convolution  In  tlna  do- 
aata  la  nultlpllcntlon  In  frequency  domain. 

hr  ns ample; 

*<*>  * «(t)  - h(t) 

ar  Kf)  C(f)  • H(l) 

•»  C(f)  - nitL 

KO 

ar  C<t>  - F'^H(f)  J 

Theoretically,  tharoforo,  deconvolution  pro- 
caduta  can  ba  uaad  to  provide  the  tranafar 
function  of  the  ayatoa. 

* more  realletlc  nodal  alallar  to  that  prasan- 
tad  In  ( 6 ) for  ultrasonic  Inspection  la  aa 
shown  to  Fig.  3. 


Into  the  values  of  the  phase  angle  profile. 
Their  study  Indicated  that  the  phase  eigna- 
tura  revealed  some  details  that  could  be  used 
to  characterise  the  flaw.  Rose.  Hast  and 
Walker  (9 ) observed  that  the  phase  profile 
elope  will  vary  for  various  flaw  cluster  situ- 
ations and,  therefore,  can  be  used  for  deter- 
mining flaw  type,  l.e.,  cluster,  porosity,  etc. 
They  also  utilised  phase  signatures  to  Identify 
disks  from  spheres  and  were  able  to  do  so  In 
all  the  cases  studied.  Rose,  Mast  and  Slklas 
(10)  also  present  a sample  problem  for  thick- 
ness measurement  utilising  phase  angle. 

Johnson  (11)  has  Indicated  that  the  phase  In- 
formation My  prove  valuable  In  Identifying  a 
particular  type  of  surface  geometry. 

Phase  signature  techniques  nay  prove  useful  In 
transducer  characterisation  and  also  In  pattern 
recognition  work  presented  In  the  next  section. 

4.  Value  of  Phase  In  Learning  Machine  Analysis 


0(t)-g(t)»p(t)*c(t)*m(t)*f(t)*m(t)*c(t)*p(t) 


Fourier  transform  of  the  above  yields 
0(f)  - P(f)  A(f ) 


A(f)-C(f)  P2(f)  M^f)  C*(f) 

F<*>  “ Off) 

A(f) 

Since  A(f)  la  In  the  denominator,  any  low 
amplitude  frequency  component  In  the  reference 
Signal  will  result  In  a considerably  noisy 
output.  The  use  of  windowing  functions  (71 
(triangular,  Hanning,  Hamming)  Is  necessary  In 
order  to  cut  out  very  low  and  very  high  fre- 
quency component.  F(f)  'a  obtained  by  point  by 
point  devlelon  of  the  two  signals  multiplied 
with  a windowing  function  and  then  Inverse 
Fourier  transformed  to  tnce  uomaln.  It  Is  im- 
portant to  point  out  that  the  correct  phase 
angle  be  used  while  performing  the  inverse 
Fourier  transform.  Tine  cornaln  representation 
so  obtained  yields  the  transfer  function 

or  fk»  reflectivity  function 

for  the  flaw. 

Implements! Ion  of  this  deconvolution  technique 
la  tedious  and  the  vaults  attained  are  not  at 
accurate  as  desired  due  to  the  nonllnearit let 
present  In  the  ultrasonic  test  aystca.  Addi- 
tion . I work  on  this  subject,  coupled  with  a 
suitable  transducer  acceptance  criteria  will 
be  carried  o t short  < to  advance  the  state  of 
the  art  on  this  subject. 

Phase  Signature  Techniques 

r v and  Schlarm  (8)  concocte.'  s theoretical 
study  on  tingle,  dual  and  trl  cylindrical 
porosity  cluster  mo  el  to  »or  no  nme  • 


Tha  flaw  sorting  study  was  carried  out  recently 
by  Rose  (12)  that  sarved  as  a feasibility  study 
In  the  sorting  of  23  selected  flaws  by  tech- 
niques of  pattern  recognition.  Overall  index 
of  performance  for  that  particular  study  turned 
out  to  bo  92Z  In  differentiating  sharp  vs. 
smooth  edge  defect.  The  type  of  flaws  con- 
sidered In  the  study  are  illustrated  for  refer- 
ence purposes  in  Table  I and  II.  Tha  flaw 
sorting  study  vas  moat  difficult  In  separating 
singular  and  multiple  smooth  type  situations. 

In  particular,  type  3,  9,  and  10  vara  difficult 
to  aoparato  from  each  other  and  In  nany  cases 
9 and  10  appeared  as  though  they  were  sharp 
edged  defects  because  of  the  superposition 
affect  occurring  fron  the  scattered  energy  at 
the  angle  bean  transducer.  Index  of  perfor- 
mance values  in  separating  the  smooth  edge 
defects  was  very  poor.  It  was.  therefore,  de- 
cided to  Incorporate  some  phase  analysis 
features  that  eight  be  useful  in  improving  the 
index  of  performance  in  this  classif icatior. 
problem.  First  attempts  were  race  at  examining 
phase  profiles  from  the  angle  beam  scattered 
data  In  two  parts,  one  fron  the  longitudinal 
component  of  the  signal  and  the  other  from 
shear  component  of  the  ultrasonic  signal. 

Data  acquisition  techniqua  is  illustrated  In 
rig.  4.  Visual  observation  of  the  phase  pro- 
file proved  tedious  and  not  immediately 
beneficial.  Here  advanced  vork  in  pattern 
recognition,  however  in  using  computerised 
feature  extrac.ron  anu  data  analysis  might 
prove  worthwhile.  It  was  decided  though  to 
uaa  normal  bean  transducer  to  get  sore  Idea  If 
the  phase  profile  could  be  used  in  this  case. 
Keeping  In  mind  that  tha  angle  beam  approach 
was  used  to  leec  at  longitudinal  and  shear 
component  coning  from  the  mode  conversion 
occurring  at  the  flaw  In  question,  phvsics 
dictates  that  a normal  beam  transducer  night 
be  suitable  for  flaw  type  determination,  in 
particular,  when  examining  the  phase  profile. 
Sample  results  of  the  normal  beam  data 
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•oqulsltlon  (ulalfM  ar*  shown  la  rig.  S. 

Flaw  Multiplicity  caa  be  caally  determined  by 
******  lag  tha  nuabar  of  dlacontlaultlaa  la  phaaa 
prof 11a  ovar  aa  actlva  frequency  raaga.  It 
night  ba  polntad  out  that  a caraful  examination 
of  slope  chan go*  la  frequency  aaplltuda  proflla 
could  ravoal  tha  Identical  information  but  ob- 
aarvatloa  of  slopa  changaa  la  nuch  acre  diffi- 
cult whan  coaparad  to  tha  largo  dlaeoatlnultlaa 
that  occur  la  tha  phaaa  proflla  In  Fig.  6. 

Additional  work  calling  for  coaplotaly  autoaatad 
computerised  feature  extraction  and  than  In- 
corporation Into  tha  various  tacholquea  In 
pattern  recognition  described  by  Rosa,  namely, 
that  of  tha  nearest  neighbor  rule  concept,  the 
Flachar  linear  dlscrlalnant  function  aad 
application  of  adaptive  learning  network  will 
certainly  be  worthwhile  If  Indeed  this  theore- 
tical approach  becomes  a practical  approach  to 
aolvlng  problem*  In  flaw  classification. 
Certainly  tha  results  presented  hare  appear 
proalslng  but  algorithms  combining  fusty  logic 
technique  with  tha  advanced  concepts  of  pattern 
recognition  are  definitely  required  along  with 
a aultable  transducer  acceptance  criteria  on 
whether  the  phase  data  should  bo  Included  In 
tha  analysla  or  not. 
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table  2-jest  Specimen  Flaw  Characteristics 
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Fig.  2 - A Basic  Deconvolution  Model 
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Fig.  3 - A Model  for  Ultrasonic  Test  System 
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Fig.  4 - Data  Acquisition  Technique 
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APPENDIX  B 

Article  1.  "Transducer  Compensation  Concepts  in  Flaw  Classification" 

Joaeph  L.  Rose  and  Michael  J.  Avioli,  presented  at  the  ASNT 
Spring  Conference,  April  2-7,  1978  at  New  Orleans,  Louisiana. 
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"Transducer  Compensation  Concepts  in  Flaw  Classification" 

by 

Joseph  L.  Rose 

Professor  of  Mechanical  Engineering 
and 

Michael  J.  Avioli 

Graduate  Student,  Mechanical  Engineering  and  Mechanics 

Abstract 

Flaw  classification  analysis  is  quite  often  strongly  influenced  by  the 
type  of  ultrasonic  waveform  that  is  generated  by  an  ultrasonic  transducer. 
One  goal  of  this  paper  is  to  Introduce  procedures  that  could  possibly  make 
flaw  classification  algorithms  become  somewhat  independent  of  certain  ultra- 
sonic waveform  characteristics  being  used  in  the  data  acquisition  procedure. 
Data  acquisition  of  ultrasonic  pulse  echo  signals  depends  quite  strongly  on 
many  teat  system  characteristics,  in  particular,  special  characteristics  of 
the  ultrasonic  transducer  and  pulser-receiver  instrument  characteristics. 

A transducer  compensation  procedure  is  presented  in  this  work  that  requires 
a suitable  reference  signal ■ containing  "noise"  contributed  only  by  system 
components  external  to  the  unknown  flaw,  and  in  software,  a processing 
scheme  is  designed  to  remove  external  effects,  therefore,  allowing  concen- 
tration on  flaw  characteristics  contained  within  the  ultrasonic  signal. 

The  processing  scheme  has  four  general  components:  Acceptor,  Compen- 
sator, Comparator,  and  Evaluator.  The  acceptor  is  basically  a gate  that 
decides  whether  or  not  a particular  transuucer  is  usable  for  the  problem 
at  hand.  The  compensator  implements  a mathematical  deconvolution  process. 
The  conparatoi  uoes  a feature  by  feature  similarity  check  on  the  desired 
t«nal  and  the  compensated  signal.  The  evaluator  is  any  scheme  that  can 
lit  tern ine  the  per fomatire  of  a given  transducer.  In  particular,  the  al- 
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gorithm  under  study  may  be  used  to  evaluate  transducer  performance.  The 
evaluation  stage  is  followed  by  an  examination  of  the  comparator  results. 
Tolerances  relating  to  the  acceptability  of  a transducer  are  obtained 
through  this  final  stage. 

Model  analysis  Is  used  to  study  the  compensation  problem.  A Layered 
Model  i.  used  with  various  level,  of  system  "noise"  being  introduced,  in 
order  to  examine  the  "noise"  effect,  in  the  deconvolution  computation 
process.  Promise  for  attaining  success  in  this  difficult  compensation 
problem  is  good,  particularly  when  considering  signal  averaging  as  a signal 
processing  tool. 
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INTRODUCTION 


The  subject  of  reliability  in  ultrasonic  examination  is  tremendously 
important.  Higher  reliability  values  have  been  associated  with  sample 
problems  in  ultrasonic  inspection,  for  example,  in  thickness  measurement 
and  in  reflector  location  analysis.  The  subject  of  flaw  classification, 
on  tne  other  hand,  usually  based  on  various  features  of  pulse  shape  have 
created  various  problems  in  classification  reliability.  Errors  in  many 
cases  have  been  attributed  to  many  different  transducer  pulse  form  varia- 
tions, coming  about  because  of  such  instrumentation  changes  as  pulser, 
receiver,  and  gain  settings.  Cable  influences  have  produced  changes  in 
pulse  shape.  In  fact,  pulse  shape  variations  have  been  reported  for  the 
same  transducer  on  different  days  and  at  different  times  because  of  both 
environmental  and  other  unknown  changes  in  the  ultrasonic  instrumentation 
network.  The  goal  of  this  paper  is  to  tackle  the  problem  of  improved 
reliability  by  way  of  establishing  a suitable  transducer  compensation 
system  to  account  for  the  changes  in  ultrasonic  pulse  shape  for  one  reason 
or  another. 

Once  guidelines  for  transducer  compensation  are  established,  work  on 
the  establishment  of  high ’ reliability  flaw  classification  algorithms  can 
fce  carried  out  as  a function  of  a single  transducer  pulse  shape.  The 
Qurcrous  variations  of  transducer  pulse  shape  need  not  be  considered  in 
the  already  difficult  algorithm  development  program  for  flaw  and  material 
variation.  Variations  in  pulse  shape  will,  therefore,  be  accounted  for  by 
employing  a transducer  compensation  routine.  Obviously,  transducer  compen- 
sation will  not  always  be  possible.  A transducer  acceptance  criteria  and 
pens.! lion  network  is,  therefore,  prop  >red  in  this  paper  that  will  advance 
V state  of  the  art  in  many  ultrasonic  inspection  programs  associated  with 
advanced  work  in  pattern  recorn! -Ion  md  computer  learning  analysis  for  flaw 
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deflector  gcoretrv  analysis  is  stron^lv  influenced  h«  t!ie  tvpe  of 
ultrasonic  waveform  that  is  generated  bv  an  ultrasonic  transducer.  One 
al  of  this  paper  is  to  introduce  procedures  tnat  could  nossiblv  rake 
flaw'  classification  algorithms  become  somewhat  Independent  of  the  ultra- 
sonic cavefoms  used  in  data  acquisition  procedures,  /coulsition  of 
ultrasonic  pulse  data  depends  strongly  on  ranv  test  s”ster  characteristics, 
in  particular,  those  of  the  ultrasonic  transducer  and  pulscr-receiver 
instrunent . 

Frederick  and  Seydel  fl]  have  discussed  a conputer-based  data  process- 
ing system  and  compensation  technique  to  increase  the  effective  bandwidth 
of  an  ultrasonic  transducer.  Their  study  included  the  use  of  a reference 
function  that  aided  in  the  preservation  of  phase  information  and  the  means 
for  a point  l>v  point  complex  division.  This  technique  improved  axial  reso- 
lution, but  was  not  anplied  to  transducer  compensation  analysis.  Mucciardi 
’2,  3]  useu  deconvolution  (transducer  compensation)  as  a preprocessing  tech- 
nique ir.  his  pattern  recognition  studies.  This  v.-orh  presents  a transducer 
compensation  procedure  chat  uses  a concept  of  a transfer  function  (4]  and 
a method  of  complex  division  that  rclv  on  suitable  reference  signals  con- 
taining abborrations  contributed  only  bv  svstem  components  external  to  the 
flaw  system.  This  procedure,  therefore,  allows  concentration  on  flaw  char- 
acteristics contained  within  the  ultrasonic  signal. 

Limited  success  in  utilization  of  this  compensation  technique  -was. 
demonstrated;  the  inverted  amplitude- time  profile  being  very  noisy.  Model 
analysis  is  used  in  this  paper  to  study  Che  compensation  problem.  This 
approach  simplifies  the  generation  of  the  data  sets  the  analyst  must  use; 
whereas  the  complete  and  varied  data  set  generation  would  be  virtually  ira- 
, • ‘;sibl  fer  the  experimentalist,  because  of  the  unavailability  of  a sig- 
ntf  'cant  number  of  flaw  specimens  required  in  the  compensation  and  pattern 
recognition  studies.  A layered  model  is  used  with  the  various  levels  of 
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system  noise  being  introduced  in  order  to  examine  the  "noise"  effects  in 
the  deconvolution  computation  process.  Promise  for  attaining  success  in 
this  difficult  compensation  problem  is  good,  particularly  when  considering 
signal  averaging  as  a signal  processing  tool. 
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Theoretical  Approach 

a)  basic  Concepts 

Transducer  compensation  is  essentially  tlie  rer.oval  of  the  effects  of 
different  transducers  from  a I’.lven  test  situation.  Consider  eke  procedure 
for  determinin'’,  the  shape  of  a subsurface  flaw  by  ultrasonic  methods. 
Physicallv,  a transducer  is  excited  by  a voltage  spike  (G(o>)  P(o))).  which 
distorts  the  crystal  in  the  transducer  (T^(u>)).  This  distortion  is  propa- 
gated as  a pressure  wave  through  the  medium  hosting  the  flaw.  Upon  incidence 
on  the  flaw,  the  wave  is  reflected  distorting  the  transducer’s  crystal 
IT,W),  which  generates  a time-varving  voltage.  This  voltage-time  profile 
is  amplified  and  is  either  displayed  on  a CUT,  digitized  and  stored,  or  both. 

The  host  medium,  the  flaw,  and  counlant  ray  he  viewed  as  a physical 
svsten  excited  by  the  initial  distortion  of  the  transducer  crvstal  and  re- 
sponding with  another  distortion  of  the  crystal.  Mien  transducers  are  changed, 
crystals  are  changed,  and  therefore  excitations  and  responses  from  the  same 
nhvsical  system  are  changed. 

The  linear  systems  approach  to  this  problem  is  to  consider  the  responses 
of  the  flaw  system  as  the  convolution  of  the  excitation  and  impulse  response 
of  tlie  system.  Mote  that  it  is  necessary  to  assure  that  on  receive  the  trans- 
ducer is  linear.  This  is  not  unreasonable  since  signal  levels  are  nviite 
small  or.  receive. 

Mathematically, 

J*f(t>  p,(r-  t)  dt^»  F(o>)  • G(aO 

- OO 

.’here  F(o»)  * Fourier  Transform  of  t(t) 

G(h>l  Fouriei  Transform  of  C.(  tl 

i i>  system  dta;rar  tor  the  nV.  w situation  appears  in  Fig.  la. 


NAKC-GSED-121 
Page  37 


IjCcu)  * F(o>)  ■ (u») 


where  I ^ (o>)  ■ Fourier  Transfom  of  input  i 

P^(w)  • Fourier  Transfom  of  output  i 

S(o»)  ■ Fourier  Transfom  of  Fouplant,  Vediun  and  Tlaw 
(known  as  the  s'’sten  transfer  function) 

This  model  adnits  a different  definition  of  transducer  compensation. 

Transducer  compensation  is  the  determination  of  the  system  transfer  function 

(which  is  considered  fixed). 

Frenuency  domain  analysis  shows  that  for  two  different  transducers, 
the  following 


S - 0. 


*2 


Implying  that  the  system  transfer  function  is  independent  of  transducer, 
h)  Deconvolut ion 

Deconvolut ion  is  the  process  fcv  which  th«  transfer  function  of  a physical 
system  is  determined.  Analytically,  this  is  a point  by  point  division  in 
the  frenuenev  donain,  0/1. 

T ie  output  0 is  usually  taken  as  the  return  echo  from  the  test  piece. 

The  input  is  determined  by  looking  at  the  return  echo  from  the  hack  vail  of 
.i  test  Mock  made  of  the  same  material  as  the  system  under  test. 

The  Fourier  Transforms  of  O and  I are  taken  and  complex  division  imple- 
mented. The  result  of  this  division  is  the  system  transfer  function. 
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Generalized  Transducer  Compensation 
Pattern  recognition  algorithms  should  be  designed  with  respect  to  a 
single  well  characterized  transducer  during  the  algorithm  development 
stage.  Practicality,  however,  requires  that  existing  algorithms  be  able 

f - 

to  handle  inputs  from  a various  assortment  of  transducers.  The  algorithm 
development  process  is  a long  and  tedious  task.  One  approach  is  to  em- 
ploy the  so-called  "shot  gun"  method.  This  method  involves  obtaining  data 
from  a multitude  of  transducers  and  designing  an  algorithm,  that  in  some 
average  sense,  performs  well  for  this  collection  of  transducers.  Unfortu- 
nately, a "typical"  set  of  transducers  does  not  exist. 

The  alternative  to  the  "shot  gun"  approach  is  to  implement  a prepro- 
cessing routine  that  would  modify  the  outputs  of  a set  of  transducers 
that  have  passed  an  acceptance  criteria  to  "look"  like  they  came  from  the 
algorithm  design  transducer.  Figure  2 shows  in  flow  chart  form,  the  gen- 
eral concept  of  compensation.  Here  it  is  assumed  that  a particular  trans- 
ducer, in  this  instance,  transducer  K,  was  used  to  develop  an  algorithm 
with  high  reliability  of  performance.  Transducers  labeled  1 through  N 
(excluding  K)  are  first  evaluated  on  a go  no-go  basis  by  an  acceptance 
criteria.  Those  surviving  this  test  are  then  used  to  probe  the  flaw  and 
material  system.  Each  transducer  will  produce  a different  output.  These 
outputs  are  entered  into  a "transducer  compensation"  routine  in  order  to 
simulate  the  output  K (input  to  the  algorithm).  Not  all  transducers  will 
be  amenable  to  this  operation.  The  body  of  this  paper  will  demonstrate 
chat  bandwidth  and  noise  characteristics  are  critical. 

I 

I 
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Acceptor  Design  Concept 

Techniques  for  designing  the  logic  associated  with  the  development  of 
an  acceptor  that  would  be  useful  in  pattern  recognition  analysis  are  out- 
lined in  this  section.  It  is  proposed  to  make  use  of  the  theoretical  de- 
convolution  routine  presented  in  this  paper  to  establish  guidelines  for 
the  acceptor  design  which  essentially  just  provides  us  with  a go  no-go 
criteria  for  a new  transducer.  Let  us  consider  a modified  diagram  of  the 
generalized  transducer  compensation  concept  that  was  presented  earlier. 

See  Figure  3.  Design  input  for  the  acceptor  will  come  about  by  examining 
the  output  functions  after  deconvolution  or  by  examining  index  of  perform- 
ance for  a set  of  data  after  passing  through  the  high  reliability  algorithm. 
In  either  case,  experience  based  on  both  theoretical  and  experimental 
analysis  will  be  used  to  establish  an  acceptance  criteria  for  a given 
transducer.  It  is  essential  in  this  approach  that  a well  characterized 
transducer  be  used  with  a list  of  characteristic  features.  Rather  than 
carry  out  a detailed  study  of  a fairly  large  number  of  transducers,  it  is 
proposed  to  consider  theoretical  variations  on  the  transducer  pulse  form 
from  which  the  high  reliability  algorithm  was  designed.  The  acceptor  de- 
sign will  depend  significantly  on  the  system  "noise".  A procedure  is 
proposed  in  this  paper  that  combines  various  samples  of  system  "noise" 
with  the  different  input  functions  in  an  attempt  to  establish  guidelines 
on  transducer  pulse  forms  that  would  obtain  reasonable  output  results  by 
way  of  the  deconvolution  routine.  Acceptor  design  can,  therefore,  be 
based  on  a theoretical  transfer  function  with  the  system  "noise"  providing 
us  with  the  major  errors  or  acceptance  limits  in  the  deconvolution  process. 
Note  that  a noise  free  system  indicates  that  all  transducers  are  acceptable 
I ’it'1'  the  deconvoluLion  process  always  performs  satisfactorily  in  this  line- 
systems  approach. 
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Theoretical  Deconvolution  Sample  Problems 
The  general  scheme  employed  in  the  following  example  is  shown  in 
Figure  4.  The  elements  labeled  "noise"  are  all  distinct  from  each  other 
and  for  this  study  had  three  possible  states:  no  noise,  one  shot  noise, 
and  noise  averaged  thirty-two  times.  The  summing  Junctions  represent  the 
fact  that  additive  noise  was  assumed  in  these  examples. 

Realistic  input  ultrasonic  waveforms  were  selected  for  this  study. 

A Layered  Media  Transfer  function  was  taken  form  a closed  fom  solution 
/ r \ 

in  Brekhovskikh v J as  outlined  in  TaLle  I and  illustrated  in  Fig.  9a. 

Two  different  ultrasonic  waveforms  were  selected  as  input  to  the  three 
layered  structure.  The  input  function  characteristics  are  summarized  in 
Table  II  and  shown  in  Figures  5a,  5b,  6a,  and  6b.  Output  functions  util- 
izing the  layered  media  transfer  function  are  presented  in  Figures  7a,  7b ■ 
8a,  and  8b.  The  transfer  function  for  these  functions  was  calculated  with 
the  deconvolution  calculation  process. 


These  functions  are  shown  in  Figures  9b  and  9c.  The  echoes  were  obtained 
by  first  calculating  TF^  and  TF?  and  then  multiplying  the  rourier  Transforms 
of  innut  one  and  input  two  by  TF^  and  TF,  respectively,  and  finally  inverting 
the  products. 

It  was  noticed  that  the  transfer  functions  obtained  bv  deconvolution 
v re  somewhat  distorted  and  shifted  to  the  left.  This  is  attributed  to 
is  Ion  Ir.orithn  noise  (see  Appendix  I).  The  waveforms  obtained  through 
.)  . i i n wcr<’  not  distorted  and  were  essentially  as  in  Figures  5c  and  6c. 
eld  seen  that  once  the  transfer  function  of  the  system  was  obtained, 
i ie  would  be  no  m.od  t calculate  the  tire  domain  echo,  ’’any  ultrasonic 
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tests,  however,  are  based  on  such  parameters  as  peah-to-peab  value,  pulse 
duration,  fall  tires  (decay)  and  so  forth.  These  values  are  not  readily 
available  fron  the  transfer  function,  although  studies  na”  he  made  to 
deduce  the  relations  between  certain  useful  parameters  and  their  manifes- 
tation in  the  transfer  function.  Therefore,  in  order  to  present  results 
in  a format  convenient  to  both  the  mathematician  and  engineer,  anpltude 
versus  time  analysis  is  covered  along  with  the  more  basic  transfer  function. 
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Deconvolution  with  Svstem  "Noise" 

System  "noise"  was  added  to  the  theoretical  results  presented  in 
Figures  5a  and  6a.  A different  set  of  "noise"  was  added  to  each  signal, 
a sample  of  which  is  shown  in  Figures  10a  and  10b.  "Noise"  was  produced 
from  fast  data  acquisition  consisting  of  UTA2  pulser-receiver  unit,  step- 
less gate,  and  Biomation  8100  A/D  converter,  as  data  was  transferrred  into 
a PDP  11/05  minicomputer.  The  goal  of  this  study  was  to  evaluate  the  "noise" 
influences  on  the  deconvolution  process.  First  of  all,  let's  consider  the 
process  of  making  the  output  function  number  2 look  like  output  function 
number  1. 


V 
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This  result  is  shown  in  Figures  11c  and  12a.  The  result  is  quite  noisy, 
but  it  is  obvious  that  the  result  resembles  the  ideal  function  0 , shown 
in  Fig.  7. 

The  reverse  computation  for  O'  is  shown  in  Figures  lib  and  12b. 


This  result  resembles  the  ideal  result  shown  in  Fig.  8. 

It  should  be  noted  that  the  result  0^'  is  substantially  "cleaner"  than 
result  O^'.  This  is  the  effect  of  bandwidth.  The  0^/1^  complex  division 
was  essentially  "narrow  band",  while  the  0^/1^  division  was  essentially 
"broad  band". 


The  inverse  of  the  power  spectrum  of  a pulse  basical  Iv  looks  like  a 
parabola.  In  the  frequency  domain,  this  means  luw  and  high  frequencies  arc 
accented.  The  broader  or  wider  the  "parabola",  the  lesser  the  effect  on  the 
i.  and  Inga  frequency  portions  of  the  transfer  function. 

Figure  13  gives  .•  pictorial  representation  of  the  complex  division 
process . Complex  numbers  may  be  repre  .ented  in  polar  form.  In  this  form, 


I 
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15c,  and  16b. 


V - 


X1  °2 


V ‘ 


X2  °1 


Results  in  this  case  are  excellent,  which  can  be  seen  by  comparing 
results  in  Figure  16a  and  16b  with  those  in  Figures  7 and  8,  respectively. 


the  spectrum  is  the  modulasNRe  + Ini  . The  other  part,  the  argument, 
tan  X (Im/Ke)  is  called  the  phase  angle.  Spectrums  are  used  because 
they  give  a view  of  complex  division  that  is  nore  easily  grasped  than 
the  use  of  the  complex  plane.  The  phase  angles  are  used  in  the  compu- 
tation process,  but  are  not  shown  here. 

Let  us  now  consider  some  improvements  in  the  results  presented  in 
Figures  lib,  11c,  12a,  and  12b.  Signal  averaging  will  be  used.  A sample 
"noise"  profile  obtained  from  an  ensemble  of  32  waveforms  is  illustrated 
in  Figures  14a  and  14b.  Results  for  the  transfer  function  and  the  cor- 
responding inverse  deconvolved  waveforms  are  shown  in  Figure  15b,  16a, 
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The  study  of  transducer  compensation  presented  in  this  paper  shows 
great  promise  for  solving  many  difficult  problems  In  ultrasonic  analysis. 
Obviously,  the  work  on  combining  theoretical  deconvolution  concepts  with 
the  "noise"  problems  associated  with  real  ultrasonic  test  systems  calls 
for  expanded  work  efforts  and  analysis.  Rattier  than  use  a layered  media 
transfer  function  as  illustrated  in  this  paper,  it  will  be  necessary  to 
obtain  real  material  and  flaw  transfer  characteristics  with  very  careful 
and  clever  ultrasonic  data  acquisition  techniques.  The  transducer  compen- 
sation program  of  study  presented  here  can  be  used  to  establish  suitable 
limits  on  various  features  associated  with  transducer  characterization,  as 
Illustrated  in  the  discussions  on  comparator  and  acceptor. 

The  linearity  assumption  in  our  systems  model  assumes  that  the  pulser 

and  receiver  units  in  a flaw  detection  instrument  can  be  operated  completely 

Independent  of  each  other.  An  evaluation  of  the  limits  and  applicability  of 

this  assumption  calls  for  a complete  understanding  of  all  controls  associated 

with  a flaw  detection  instrument.  The  quality  of  the  noise  considered  in 

this  study  is  of  considerable  importance,  the  "noise"  in  our  case  being 

associated  with  our  particular  fast  data  acquisition  system  utilizing  a 

• « 

ITTA2  pulser-receiver  and  a Biomation  8100  analog  to  digital  converter.  An 
evaluation  of  transducer  compensation  performance  for  other  test  systems 
and  noise  components,  however,  can  certainly  be  carried  out  with  the  pro- 
gram of  study  presented  in  this  paper. 
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Comutcr  Deconvolution  Procedure 


Consider  the  complex  multiplication 


(a  + jb)(c  + jd)  - (e  + jf) 


.-here  a.  b,  e,  and  ? arc  knm.*n  values , c and  d are  to  he  found  (transfer  function) 


ac  - bd  ■ e 


(real  parts  enual) 


he  - ad  • f 


(imaginary  parts  enual) 


Rewrite  as 


[:  1 [:]•[;] 


2 2 

then  if  a + b ♦ 0 (Determinant  non-zero),  the  solution  for  the  system  is 


*2-2  [-*  J U'[J 


(uniouely) 


The  literals  nay  be  considered  as  subscripted  variables  or  array  variables. 


The  values  for  c and  d mav  be  calculated  each  time  a”  + b"  ■ 0.  The  senutnees 


2 2 

c and  d now  are  determined  except  for  "holes"  at  the  locations  where  a"  + b-  + 0.' 


A mathematical  argument*  requires  that  three  "holes"  he  filled  with  interpolated 


values.  The  interpolation  schene  then  becomes  critical.  It  seems  the  only 


available  criteria  to  date  is  one  based  on  comparisons  letween  analytically 


venerated  transfer  functions  and  ones  obtained  with  Che  above  process  using 


a particular  interpolating  scheme. 


Method  1 


Method  2 


I Method  3 


! i‘  e = a,  f " b,  the  division  of  two  identical  complex  numbers,  although 


i » ’>’  » 0 at  tines  would  require  the  continuous  result  1. 
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Tania  I — rireKhovsI'.lkh  Layered  Modi  a Transfer  Function  Paraneters 


density 

Sr./cm^ 

v.-ave  velocity 
cn/  sec 

thickness 

cm 

layer 

1 

2.71 

.63 

.^336 

layer 

2 

1.18 

.27 

.0168 

— 

layer  3 2.71 

.63 

.0336 

fj 

l 

j£j».l«_  11  - Sariple  Ultrasonic  Input  Vavefoms 

Obtained  from  Real 

Experimental  Data 

■: 

|0 

6 d)J  down 

6 dR  down 

Center  Frequency 

MHz 

. . _ 

Kandvidth 

Center  Frenuencv 

1 

MHz 

MI!z 

i 

Function  ;?1  3.81 

' 

3.22 

3.47 

s 

] 

.1 

Function  *2  5.18 

2.05 

5.32 

Fig.  la  - A Basic  Deconvolution  Model 
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Fig.  12a  - Deconvolved  Result  with  Noise 

(making  Output  Function  Two  look  like  Output  Function  One) 


r l-'l  - Deem  'olved  Result  with  Noise 

(‘"’'•‘"8  °utP"1  Function  One  look  like  Oi  tput  Function  Two) 
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Low  Frequency 
D.C.  Levels 

Noise 


High  Frequency 
Noise 


Distorted  Mid-Kange 


Fig.  13  - Illustration  of  the  Deconvolution  Process 
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Time  {usee) 


Fig.  14a  - Sample  of  Averaged  Noise  1 


Time'^/sec) 


Fig.  14b  - Sample  of  Averaged  Noise  2 
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"if  1 bb 


'■'•crnvol'  <i  Result  with  Averaged  Noise 

("1  tug  Oj.put  Function  One  look  like  Output  Function  Two) 


I I 

Frequency  (>tiz) 


' - . 5;  - lip  i run  of  Input  One 


Frequency  (MHz) 

9a  - Theoretical  Transfer  Function  for  a Sample  Layered  Media  Problem 


Fig.  9b  - 


I 1 " 1 

Frequency  (Jfliz) 


Deconvolved  Transfer  Function  (O./I.) 
(shouinc  mathematical  noise) 


A 

A / \ 


Fie. 


! ! ! <iii 

• i i i , i . 

Frequency  (Mllz)  IQ 

- Deconvolved  Transfer  Function 
(showing  mathematical  noise) 


'Ill 


I 


i r 


Time  (/Lisec) 

Fig.  10b  - Sample  Noise  2 


Time  (fisec) 

Fig.  12a  - Deconvolved  Result  with  Noise 


(making  Output  Function  Two  look  like  Output  Function  One) 


Time  fyxsec) 

Fig.  12b  - Deconvolution  Result  with  Noise 

(making  Output  function  One  look  like  Output  Function  Two) 


Tirce  (/nsec) 

Fig.  14b-  Sample  of  Averaged  Noise  2 


Fig.  15b  - Deconvolved  Transfer  Function  (CL /I,) 


Fig.  i;>c  iV-  ".tolveii  Transfer  ,'uretion  (0_.  1 ) 
(snowing  averaging  effect  on  noise)1 


Time  (jLisec) 


Fig.  16a  - Deconvolved  Result  with  Averaged  Noise 

(caking  Output  Function  Two  look  like  Output  Function  One) 


Tima  (/asec) 


'.fib  - Deconvolved  Result  with  Averaged  Noise 

< nukinc  Output  Function  0„e  look  like  Ounut  Function  Two) 

I! 


